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Abstract. Tycho's supernova remnant (SNR) was observed during 1997 and 1998 with the HEGRA Gerenkov 
Telescope System in a search for gamma-ray emission at energies above ~ 1 TeV. An analysis of these data, ~65 
>v> " hours in total, resulted in no evidence for TeV gamma-ray emission. The 3a upper limit to the gamma-ray flux 

%^ , (>1 TeV) from Tycho is estimated at 5.78 x lO"'^^ photons cm~'^ s~^, or 33 milli-Grab. We interpret our upper 

j^ ' limit within the framework of the following scenarios: (1) that the observed hard X-ray tail is due to synchrotron 

emission. A lower limit on the magnetic field within Tycho may be estimated B > 22 nG, assuming that the 
RXTE-detected X-rays were due to synchrotron emission. However, using results from a detailed model of the 
ASCA emission, a more conservative lower limit B > 6 /iG is derived. (2) the hadronic model of Drury, Aharonian 
& Volk, and (3) the more recent time-dependent kinetic theory of Berezhko & Volk. Our upper limit lies within 
the range of predicted values of both hadronic models, according to uncertainties in physical parameters of Tycho, 
and shock acceleration details. In the latter case, the model was scaled to suit the parameters of Tycho and re- 
normalised to account for a simplification of the original model. We find that we cannot rule out Tycho as a 
potential contributor at an average level to the Galactic cosmic-ray fiux. 

Key words. Gamma rays: observations - ISM: supernova remnants: individual objects: Tycho's SNR 

1. Introduction of 7-rays in SNRs is thought to result from the interac- 
tion of shock-accelerated particles (hadrons and electrons) 

The search for gamma-ray emission of TeV energies from ^^^^ ^^^ interstellar medium (ISM) and local soft photon 

supernova remnants (SNRs) m recent years is motivated g^j^^^ ^^^^^ relatively unattenuated over long distances 

by the need to explain the origin of Galactic cosmic-rays ^^^ preserving the production site directionality, gamma 

(CR). SNRs are long-believed primarily responsible for the radiation is one of the most accessible tracers of CR ac- 

Galactic CR popula tion, m atchin g the C R energetics and celeration in the universe, 
spectral index (Volk 1997, Baring 2000a). The production 



TeV gamma radiation is primarily expected from 
Send offprint requests to: G.P. Rowell, e-mail: two channels (1) Collisions of hadronic CRs, produc- 
Gavin.Rowell@mpi-hd.mpg.de ing gamma-rays via 7r° decay and (2) CR electrons up- 
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scattering soft photons via the inverse Compton (IC) 
process, and CR electron cohisions via Bremsstrahlung. 
Detailed modelling of SNR environments has revealed dis- 
tinct spectral features in the GeV/TeV regime for these 
processes, and in combination with radio and X-ray obser- 
vations, those in the TeV regime are deemed vital in estab- 
lishing SNRs as sites of CR production (Drury, Aharonian 
& V olk et al. |1994 Nait o fc T akahara |1994 Berezhko & 
Volk |1997| , Baring et al. |2000bD . 

So far, evidence for TeV emission has come from 
CANGAROO observations of two southern hemisphere 
SNRs, SN 1 006, and SNR R X J17 13. 7-3946 (Tanimori 

2000|) 



et al. 1998 



and the HEGRA 



Muraishi et al. 
Cerenkov Telescope CT-System (High Energy Gamma 
Ray Astronomy Cerenkov Telesco pe) a fter deep observa- 
tions of Cas-A (Piihlhofer et al. |l999 , Aharonian et al. 
2000a). Upper limits are reported for those other promis- 
ing SNR candidates observ ed so far in both t he Te V and 
PeV regim es (Buckley et al. 1998 , Prosch et al. 1996 , Goret 
et al. |l99g| , Allen et al. |l99§Rowell et al. |2000[ ), including 
Tycho's SNR, for which the Whipple Collaboration ob- 
tained 14.5 hours data. The CANGAROO and HEGRA 
detections might be interpreted in the framework of the 
SNR as a source of multi-TeV CR electrons by virtue of 
a strong non-thermal tail in their X-ray spectra above 1 
keV (Koyama et al. |l995i |1997| , Allen et al. |1995D. O n the 
other hand as argued by Aharonian & Atoyan (1999) and 



Berezkho, Ksenofontov fc Petukhov (1999) for SN 1006, 



and by Atoyan et al. (|2000|) for Cas-A, the TeV results 
do not rule out the hadron channel. It is nevertheless less 
clear at the moment as to the location of CR hadron ac- 
celerators in our galaxy. 

Tycho's SNR (G120.1-H1.4, 3C 10) is one of the most 
intensely studied SNRs. It is an archetypal shell- type 
(radio & X-Ray) SNR, formed most likely from a type 
la supernova (SN 1572), and has expan ded a t ^0.1% 



yr -'^ to a radius ^ 4' (Katz-Stone et al. 2000 and ref- 



erences therein). One estimate of distance is put at 2.2 
kpc (Albinson et al. 1986 and references therein) based 
on proper motion studies, absorption against field stars 
and the fact that Tycho appears embedded in the Perseus 
arm of the Milky Wa y. A h igher distance estimate of 4.5 
kpc (Schwarz et al. 1995) is derived from a model of 
the HI spectrum and number of HI absorbing features in 
the region. The radio synchrotron emission (20 & 90 cm) 
shows variation in photon index that appears correlated 
with edge filaments, perhaps tracing regions of enhanced 
magnetic field and particle acceleration (Katz-Stone et 
al. |200C ). Studies at the HI 21 cm line indicate expan- 
sion on the eastern side of Tycho is slowed by a region of 
higher density (160-325 cm"'^, Reynoso et al. 1999), sug- 
gesting that Tycho's SNR may not be expanding into such 
a homogeneous region as earlier believed. The X-ray emis- 
sion shows very strong line features and overall, is well fit 
by thermal brem sstrahlung components (Fink et al. 1994, 
Petre et al. 1999 ). A power law however is necessary to fit 
a hard X-ray tail. This presumably non-thermal tail above 
1 keV suggests that Tycho may a source of electrons up to 



1997 



1998 



Total 



Obs. time (hrs) 
Selected runs 
Configuration 
Zenith angle range 



19.7 44.9 64.6 

50 115 165(82%) 

CT3-6 CT3-6 

32-41° 32-45° 



CR Event Rate (Hz) 8.4 (10.3) 



8.5 



8.5 



Table 1. Summary of Tycho Observations by the HEGRA 
CT-System. The CR event rate is after the size cut of 40 
photoelectrons and the CR rate in brackets for 1997 is 
prior to a cut of 10 photoelectrons on second brightest 
image pixel. 



^100 TeV, thereby joining the growing number of SNRs 



that exhibit this feature (see Petre et al. 1999 for a re- 
view). 

2. Observations 

The observations of Tycho's SNR were made using the 
HEGRA CT-System incorporating, at the time of data 
taking (1997 & 1998), four identical imaging air Cerenkov 
telescopes|j(CT3-6). The telescopes are operated in coinci- 
dence to achieve a stereoscopic view of Cerenkov light air 
showers induced by 7 and CR primaries. The CT-System 
is situated on the Roque de los Muchachos at La Palma 
(2200 m asl, 28° 45' N 17° 54' W) and at the zenith angles 
of these data, operates at an energy threshold for 7-ray 
primaries of '--^l TeV (see Konopelko et al. 1999 and ref- 
erences therein). 

The Tycho data were accumulated in Jul-Sep 1997 
and Aug-Dec 1998 in 20—30 minute runs utilising the 
wobble mode of observation. In this mode, the source is 
offset ±0.5° in declination in alternating runs, enabling 
the background or OFF source data to be estimated from 
the opposing position (or a series of positions) 1.0° away. 
Table H gives a summary of these observations. 

Approximately 20% of all runs were rejected, due pri- 
marily to weather effects. Bad runs were identified ac- 
cording to the criteria that the image width from each 
triggered telescope and stereo trigger rate match their ex- 
pected values accordingly, within 6% for width and 15% 
for rate (representing levels < 4(7 outside their Gaussian 
distribution means). The bad runs actually comprised the 
non-Gaussian tails of these distributions. The trigger rate 
for CRs was significantly lower in 1998 compared to 1997 
(by ^--^25%), attributed mainly to photomultiplier tube 
(PMT) fatigue and mirror degradation. We accounted for 
this threshold difference by artifically raising the trigger 
threshold of 1997 data to 1998 levels by applying a 10 
photoelectron threshold to the second-brightest pixel in 



^ Since 1999, the full five telescopes have been in operation, 
CT2 being the most recent telescope brought online. 
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the image, mimicking the hardware trigger threshold, at 
8 photoelectrons. The effect of this software trigger is to 
simphfy analysis and in particular, the estimation of up- 
per limits. It may be justified given the negligible effect 
on the 7-ray threshold of this cut (following the 40 photo- 
electron cut described below) and the fact that the 1997 
observation time is ~45% that of 1998. The 7-ray energy 
thresholdn for both years' data were found from simula- 
tions to be consistent within systematic errors (15%) to a 
value of 1 TeV. 

Event filtering follows the standard method in that 
images containing more than 40 photoelectrons (size) are 
selected for further analysis. The rejection of the CR back- 
ground is achieved with cuts on the image shape and the 
location of it's reconstructed source position. The recon- 
struction is described as "Algorithm 1" by Hofmann et 
al. (1999) and the shape cut used is the so-called mean 
scaled width (msw). The msw is the average of the scaled 
widths from each accepted image. Each image is scaled 
according to it's expected value (7-ray hypothesis), de- 
pendent upon image size, zenith angle of observation and 
impact distance of the event from the CT-System centre. 
Following cuts on msw and 6, the difference between the 
reconstructed and actual source positions, the stereoscopic 
method permits unprecedented angular and energy reso- 
lution for 7-ray images on an event-by-ev ent bas is of ~ 6' 
and ^20% respectively (Aharonian et al. 2000b|) . 

We use a combination of "tight cuts" in the shape 
(0.4<77isw<l.l) and direction (6*^ < 0.02 deg^) criteria 
to preferentially select 7— ray images against those from 
CR. Tight cuts are optimal in the search for weak sig- 
nals in background dominated statistics (3— Sir level), with 
the caveat that systematic effects may become more im- 
portant due to the cuts lying close to distribution tails. 
Tycho is a marginally extended source for the HEGRA 
CT-System, and thus the optimal 0^ cut will be slightly 
larger than that for a point source at ~0.015 deg^. A sim- 
ple Monte Carlo, matching roughly the statistics of the 
background of our Tycho dataset, was used to determine 
an optimal 9^ cut, assuming various morphologies for the 
SNR with the results displayed in Fig. |l|. Given the uncer- 
tainty in the SNR morphology (disk, Gaussian, annular or 
a combination), a 6^ cut of 0.02 deg^ was deemed suitable. 

A useful indicator of data quality is the CR acceptance 
('^cr) of the m,sw cut on a run-by-run basis (Fig. 0). We 
can see that kcr is consistent over the entire dataset and 
moreover, the cut of 10 p.e on the second brightest pixel 
applied to 1997 data introduced no significant change in 

KCR- 

3. Results 

We evaluate the statistical signific ance of the post-cut ex- 
cess using equation 9 of Li & Ma (1982). The background 
or OFF source counts are taken from three control re- 
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Fig. 1. Signal to noise ratios (unnormalised curves) from 
a Monte Carlo for sources of various morphologies. The 
HEGRA CT-System point spread function (PSF, with a 
specified in the plot) is assumed to be a Gaussian and 
based on measurements on the Crab (Aharonian et al. 
pOOOd ). Here 100000 background events were sampled, re- 
producing roughly the statistics of the Tycho dataset, and 
the number of source events chosen to reproduce a ~ 5cr 
maximum excess. 
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0.02 <i^ 



-> -e 
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Run number [ arbitrary J 



The maximum differential trigger rate, for a number of pho- 
ton indices. 



Fig. 2. Cosmic ray acceptance (kcr) of the ttisw cut on a 
run-by-run basis. The 1997 data were also subjected to a 
cut on the second-brightest pixel of 10 p.e. 



gions at position angles 180° & ±90°, on a circle 0.5° 
radius centred on the tracking position, as the wobble 
mode of observation permits. A normalisation factor of 
i is therefore used in the significance calculation. By us- 
ing more than one background region, effects of skynoise, 
camera response and zenith angle difference are averaged 
out further (on top of that already afforded by the stereo- 
scopic technique) , reducing somewhat any systematics ex- 
pected from such properties, and moreover, reducing the 
statistical error in the background. A distribution of (P 
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Fig. 3. Distribution of 6^ for combined 1997 and 1998 ob- 
servations of Tycho's SNR (ON source, solid line). The 
OFF source data (dashed line) are taken from three posi- 
tions in the field of view (see text) and scaled by a factor 
i . The optimal cut is indicated by the dotted vertical line 





1997 








1998 


Total 


ON 


127 








259 


386 


OFF (EffeO 


383 








871 


1254 


ON-OFF (a) 


-0.7 








-1.6 


-1.4 


3a UL (Crab) 












0.033 


3a UL (F(> ITeV) x 


10" 


-13 ph 


cm 


-s-i) 


5.78 



Table 2. Post-cut (after msw & 9^ cuts) statistics of the 
Tycho data (significance calculated using equation 9 of 
Li & Ma (1983)) taking as the background, the sum of 
positions 6^. The 3a upper limit (99.7% UL) is expressed 
in Crab and absolute units (see text). 



for the combined data after application of the msw cut 
discussed above is presented in Fig. y. We summarise 
further the results in Table g. No significant excess was 
seen in either dataset, and the combined excess repre- 
sents a significance of — 1.4ct. We may express our upper 
limit in absolute units by using the TcV flux from the 
Crab Nebula, thereby taking advantage of the substan- 
tial simulation effort in deriving the absolute Crab flux 
(Aharonian et al. |2000bD Flu ctuating the ON-OFF ex- 
cess according to Helene (1983) and using the method de- 
scribed by Aharonian et al. |20004 our 3cr UL (99.7%) 
corresponds to 33 milli-CrabPTHere, those Crab data 
used to estimate the UL (~6 hours in total), were sub- 
jected to an identical analysis as described above, and se- 
lected to match the zenith angle distribution as that for 
Tycho. Above the threshold energy of 1 TeV, our UL is 
F{> ITeV) = 5.78 x IQ-^^ ph cm^^g-i. 

Finally, we should point out that the method of 
Aharonian et al. 2000c (which downward fluctuates the 



1 Crab; F{> ETeV) = L 75 x 10'" (pf^) 



cm ^s 1 (Aharonian et al. [20001 ) 



ph 



Crab excess) to estimate the UL in Crab units was tested 
for accuracy with a simple Monte Carlo and found to 
overestimate a true 3cr UL (for the numbers of Table 
by about 25%. This overestimate is effectively canceled 
however by the negative bias of our —lAa excess, an un- 
derestimate of ~30% relative to a zero excess result. We 
therefore may assume that our 3cr UL represents closely a 
true estimate. Secondly, any integral UL or flux estimate 
quoted in absolute units will naturally depend somewhat 
on the difference in photon index between the Crab and 
Tycho. This uncertainty amounts to <30% on the UL es- 
timation for a 20% uncertainty in the values of 1 TeV and 
— 1.59 respectively for energy and photon index. 

4. Inverse Compton Interpretation 



The evidence for a non-thermal tail (Fink et al. 1994, 
GINGA data; Petre et al. |1999| , RXTE data) may sug- 
gest that Tycho's SNR is an accelerator of CR electrons to 
multi-TeV energies. In such a framework, the synchrotron 
process is assumed to account for these X-ray photons, and 
the possibility of inverse Compton (IC) TeV 7-ray emis- 
sion should be considered. However, one must be careful 
here since a power law fit is reasonable only at energies 
above ~10 keV and, alternative physical explanations to 
the synchrotron scenario do exist (see for e.g. Asvarov et 
al. [1990| , Laming |l998| , and Tatischeff et al. [19981 ). 

Assuming that the synchrotron model is valid for 
Tycho, we may use the direct relationship between the 
expected IC (fj) and X-ray energy fluxesQ] (/x) arising 
from the same electrons: 



/x(£keV) 



10 



B 



10^5 G 



(1) 



to establish a condition on the magnetic field B of the 
SNR. The mathematical caveat here is that the correct 
energy range in the TeV and keV regimes must be ad- 
hered to and that we assume that the emission regions of 
both components have the same size. In the J-function ap- 
proximation of the synchrotron emissivity for an electron, 
the IC (-B / 1 TeV) and X-ray (e / 1 keV) synchrotron en- 
ergies in Eq. (|^) are coupled according to (see Aharonian 
et aL|l997|): 



0.07 



E 



B 



1 TeV / V 10-5 G 



keV 



(2) 



so that for IC fluxes at ~1 TeV, a comparison with the X- 
ray flux at energies ~0.1 to 1.0 keV is required for values of 
B - 10 to 100 ^G expected in a SNR. Eq. (^ follows from 
the fact that the ratio of the synchrotron and IC energy 
fluxes is proportional to the ratio of the energy densities in 
the magnetic and CMB fields respectively. The dominant 
IC flux arises from the up-scattering of CMB photons by 
CR electrons of energy 10 to 1 00 Te V. Preliminary results 
from RXTE data (Petre et al. 1999) indicate a normalisa- 
tion (at 1 keV) and photon index respectively of 0.30±0.02 

" Energy flux: E^F{E) 
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photons cm ^ s ^ and — 3.18±0.02, valid at energies from 
10 to 20 keV. From earlier GINGA data Fink et al. ( |1994D 
derived a power law fit (in combination with a thermal 
bremsstrahlung and Fe line component) of photon index 
—2.72 (see Fink for error details) and normalisation of 
7.4^74 X 10^^ ph cm^^s^^ (at 1 keV), which is vahd from 
4.5 to 20 keV. A direct extrapolation to energies required 
for the B field estimate, namely 0.1 to 1.0 keV, will lead to 
overestimation of the synchrotron fiux since one would ex- 
pect a turnover near this range. A more reasonable method 
to estimate fluxes at such lower energies is to flt to the 
X-ray flux a power law spectrum with exponential cutoff. 
The RXTE data is preferred for such a fit due to the much 
reduced errors compared to those of GINGA data. We fit 
the spectrum using the function: 



E^ F{E) = AE'^-"' exp[-(S/^m) ] ergcm"^ s 



(3) 



where a (radio photon index) and E-^-i (cutoff energy) are 
free parameters and A is the normalisation determined 
from the flux density at 1 GHz S'ighz=56 Jy (Green 
2000). Eq. (P) is derived for an electron spectrum with 
exponential cutoff using the (5-function approximation for 
the electron synchrotron emissivity and is shown to ad- 
equately describe the X -ray spectrum out to energies 
~ 10^;^ (Reynolds |l999|) . Our fit yields values of a=0.7 
and i?ni=1.6 keV respectively with a differing somewhat 
from the radio photon in dex, li sted as 0.61 in Green's cat- 
alogue. Katz-Stoneet al. (2000) summarise historical mea- 
surements of the radio index as varying between 0.5 and 
0.7 at different radio bands. From their VLA study, Katz- 
Stone indicate an average index of the filaments at 0.52. 
Results are plotted in Fig. ^ showing our best fi t, and 
that when fixing a=0.61. We note that Hwang et al. (1998) 
have placed an upper limit to the non-thermal luminosity 
in the 0.5 to 10 keV regime on the order of ~10'^'* erg s^^ 
(0.5 to 10 keV), based on careful modelling of the non- 
thermal ASCA detection. In order to accommodate their 
upper limit to the non-thermal luminosity (for a distance 
of 2.3 kpc as used by Hwang) our synchrotron spectrum 
fit to RXTE data must be scaled down by a factor ~15.5, 
assuming that their upper limit is exactly lO"^"' erg s~^ 
(dotted line in Fig. ||a). 

If we take the estimated synchrotron fluxes from 0.1 to 
1 keV from the direct RXTE fit with our upper limit at 1 
TeV for the IC flux, a range of B_^ values 22 to 24 /iG, 
plotted in Fig. ^ (solid line), are obtained, which are rea- 
sonably consistent with that expected after amplification 
by the SNR shock. Using however th e lowe r synchrotron 
energetics implied by Hwang et al. ( 199^ ) we arrive at 
more conservative lower limits on _B_5 at ^6 ^G, (dotted 
line) . Such a value would result if little or no amplification 
was present. Allowing for uncertainty of one order of mag- 
nitude in the Hwang estimate (ie. 5x 10"'^'^ to 4x 10"'^'* erg 
s""'^) would, for the upper bound case, give a B field lower 
limit a factor \/4 higher, ie. ~12 /xG. 
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Fig. 4. (a) Optimal fit to the RXTE spectrum and radio 
normalisation using Eq. (^) (thin solid line) . (dashed line) 
A fit when fixing q;=0.61. (thick solid line) For compara- 
tive purposes a model fit to the XMM MOS 1 spectrum 
from Decourchelle et al. ( ^OOl ) fig 1 is included. Since in- 
terstellar absorption is not removed in this fit, the XMM 
flux below l.O keV is therefore a lower limit, (dotted line) 
Re-scaled RXTE fit to accommodate an upper limit to the 
non-thermal flux of Hwang et al. ( 199§| ). (b) (Solid line) 
Calculation of B_^ (10 /iG) of Eq. (Ill) taking as /x the 
fit given by the solid line of (a), the RXTE fit, and (dot- 
dashed line) J5_5 values obtained when using the scaled- 
down fit to a synchrotron spectrum given by the dotted 
line of (a) as suggested by Hwang et al. ( |1998 ). 



5. Comparisons with DAV 

We will consider here the production of 7-radiation by the 
hadronic or 7r°-decay channel described in two-fluid model 
by Drury, Aharonian & Volk ( |1994[ ), hereafter DAV, that 
has been used extensively in the past (see also Naito & 
Takahara 1994). The DAV prediction may be scaled sim- 
ply according to the energy budget of the SNR Esm den- 
sity of the upstream target matter n, and distance to the 
remnant d: 



1- > 1051 erg 



1 cm 



Ikpc 



(4) 
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Fig. 5. Upper limits from the present work (H-CT), and 
previous measurements (W - Whipple, H-A - HEGRA 
AIROBICC, CM - CASA-MIA. See text for references) 
with predictions of the 7r°-decay 7-ray flux from the DAV 
model. These DAV predictions from Volk 1997, use two 



limiting choices of parameters for 0, and the spectral in- 
dex (indicated on the plot). Values for the other param- 
eters were ^sn = 2 x 10^° erg ( Smith et al. |1988D , n=1.0 



cm ^ and (i=2.3 kpc (Heavens 1984). Note that a reduc- 



tion of the 'nominal theory' curve with a value 9 = 0.05 
may be argued in on the grounds that Tycho is pre-Sedov 
(see text). 



The parameter Q represents the fraction of energy avail- 
able for the acceleration of CR that are ultimately re- 
leased into the ISM. The mean value of G per SNR in 
the gal axy is estimated in the range 0.05 - 0.3 (Berezhko 
& Volk 2000), and pertains to an evolutionary state well 
into the Sedov phase when accelerated particles start to 
be released into the ISM. We shall adopt the reasonable, 
and often used value of 8 = 0.1, particle spectral index 
— 1.1 (nominal theory), and for purposes of argument that 
Tycho has progressed well into the Sedov phase. 

Fig. 1^ highlights a comparison of this DAV prediction 
under limiti ng parameter-space selections as calculated by 
Volk ( |1997[ ) with our present UL, and those of previous 
measurements (Whipple, Buckley et al. 1998 & HEGRA 
AIRO BICC , Prahl et al. |1997| , and CASA-MIA, Borione 
et al. [1995D . 

It is clear that with a factor ~4 reduction between 



the Whipple and present upper limits (assuming a spec- 
tral index of —1.1 for this comparison), the DAV model 
is now constrained when allowing for a reasonable range 
of input parameters. We are appro aching the conservative 
boundary assumed by Volk ( 1997 ) for the Sedov phase. 



An important parameter in view of the expected TeV 
emission from a SNR is it's evolutionary phase. Generally 
accepted is the notion that the maximum TeV 7-ray emis- 
sion occurs at the beginning of the Sedov phase and a 
broad maximum in the total CR energy is reached during 
the Sedov phase. At the former time the mass of swept-up 
material will just exceed that ejected, and the outer SNR 
shock begins to expand with radius R ~ t^/^ . Typical ages 
of Galactic SNRs for the onset of the Sedov phase are 500 
to 1000 yrs. Tycho is a rather young SNR (428 yrs old) ex- 
panding into an ISM with density typical of that expected 
within the Galactic plane and so it may not be clear as to 
exactly what phase Tycho is presently in. Observational 
evidence indicates that globally, Tycho's SNR is near to 
the S edov phase [R ^ t" , v = 0.46 ± 0.02, Tan & GuU 
1985 ), and is expanding into, on a global scale at least, a 
homogeneous ISM. Reynoso et al. (1997) report a similar 
value, from VLA data for the global expansion, but note 
that V varies between 0.2 and 0.8 around the shell and 
suggest the presence of denser material along the eastern 
s ide (s ee also recent 21cm observations by Reynoso et al. 
( 1999| )). The expansion rate at X-ray energies is signifi- 
cantly higher than at radio energies, but as in the radio, 
it also varies considerably azimuthally (Hughes 200C). In 
contrast to the radio behaviour however, the expansion 
rate at X-ray energies varies radially, supporting the idea 
of an ejecta-dominated, or pre-Sedov evolutionary state. 
Taking these results together with the Rayleigh- Taylor in- 
stabilities along the eastern side suggested by Velazquez 
et al. (1998), we may conclude that the global evolution 
state of Tycho's SNR is pre-Sedov, although at regions of 
high density ISM, the evolutionary state will be locally 
advanced. 

A comparison with DAV predictions should therefore 
allow for the current evolutionary state of Tycho. The 
lower curve of Fig. |g (nominal theory) may be lowered 
further since the age of Tycho is likely less than the sweep- 
up time (io=555 yrs) which signals the onset of the Sedov 
phase. A reduction in the 7-ray emissivity for pre-Sedov 
epochs can be manifested as a reduction in the 6 param- 
eter to values <0.1. A quantitative estimate for in this 
case is not trivial but generally one can expect a reduc- 
tion by a factor ~2 based on consideration of results in 
figs. 1 to 3 of DAV, describing different rates of luminos- 
ity increase with time. In the next section comparisons to 
a model dealing with time-dependence in detail are made. 



6. Comparisons with Kinetic Theory 

The most recent nonlinear kinetic models are the time- 
dependent, spherically symmetric solution of Berezhko & 
Volk 1997, hereafter BV, and the Monte Carlo simula- 



tion of a quasistationary outer SNR s hock in plane par- 
allel geometry by Baring et al. 1999 . The BV model is 
based on the nume rical solution of Berezhko, Yelshin & 
Ksenofontov ( |l994| ) and invokes a distrib ution of ejecta 



velocities (e.g. Chevalier & Liang |1989| and references 
therein) that contains very high speed components com- 
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Case 


a 


b 


c 


d+ 


BifiG) 


5 
10"^ 


30 
IQ-^ 


5 
10-^ 


5 
10-^ 


^-13 


91.00 


49.20 


20.00 


9.00 


Oe 


0.60 


0.55 


0.50 


0.50 




0.17 
15.50 


0.18 
9.00 


0.20 
4.00 


0.20 
1.80 


Qrit = 428) 


2.60 
0.06 


1.50 
0.04 


0.69 
0.04 


0.31 
0.04 



t Using a mean value for ejeeta velocities. 

Table 3. Theoretical 7-ray fluxes, FZi^ (10^-^^ ph cm^^ 
s^^) and released relative CR energies Qs from the BV 
model scaled according to parameters of Tycho's SNR (see 
text). Each case represents a choice of injection rate 77 
and B-field (labeled in the caption of Fig. ^. The pa- 
rameters subscripted with r refer to those additionally re- 
normalised to account for the effect of assuming in the BV 
model, a shock normal parallel to the B field (see results in 
Fig.g). The re-normalisation factor here is (0 = 0.1)/9s, 
and 0r(i = 428) is the relative CR energy calculated at 
the current age of Tycho's SNR. 



pared to the mean ejeeta velocity. At early evolutionary 
phases this leads to much higher shock speeds than im- 
plied by the mean ejeeta velocity, and thus to much more 
intensive CR and 7-ray production. We will in this paper 
make use of BV as it was calculated in their paper, scaling 
the parameters to those of Tycho's SNR. As we shall note 
shortly, in addition, a physical re-normalisation of these 
results is necessary. For a more complete treatment, we 
refer to a companion paper elsewhere in which the BV 
model will be calculated assuming Tycho parameters. 

A scaling of the original BV calculations is required 
since the start-up parameters (in particular ejected mass 
Mej = IOM0 and SNR energy Esn^lQ^^ erg) may differ 
from those one might expect for Tycho's SNR. For this 
scaling we may reasonably assume, that the ejeeta mass 
only effects the initial normalisation in the form of the 
sweep- up time to, and that other parameters will scale 
the flux in the same way they do the DAV prediction, ie. 
independent of time. The first row of Table Represents a nu- 
merical comparison of our upper limit with the predictions 
of the BV model, i^2i3 (^"^ units of lO^^'^ ph cm~^ s~^), 
scaled to Tycho parameters. 9s is the released fraction of 
energy available for CRs and 9 = 0.1 is the value assumed 
for an average SNR, discussed previously. Inspection of 
i^2i3 fo^' ^1^ cases in Table ^ reveals that they are incon- 
sistent with our upper limit. 

However, we must now invoke a further, physical re- 
normalisation of the Tycho-scaled BV predictions for the 
following reason. The BV model assumes a parallel shock 
geometry (shock normal parallel to the average up-stream 
B field) over the entire spherical shock surface. For spher- 
ical geometry expected of a SNR, such an assumption 
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Fig. 6. Time dependence of the 7r° decay TeV 7-ray in- 
tegral flux FZi3 (10^^"^ cm^^ s^^) from Tycho for an in- 
jection rate rj = 10~^ and B = 5/iG (dashed line, case a), 
77 = 10~^ and B = 30/xG (solid line, case b), ry = 10~* 
and B — 5fiG (dash-dotted line, case c). The dotted line 
(case d) represents a less realistic single velocity ejeeta 
case with 77 = 10~^ and B = 5/iG, as is adopted in the 
DAV model, to is the time of Sedov phase onset. These re- 
sults have been scaled according to parameters for Tycho 
(Mej = IAMq, Ssn = 2 X 10^° erg, n = 1.0 cm'^, d=2.3 
kpc. to is estimated at 555 yrs, somewhat larger than the 
current age of 428 yrs), and re- normalised by the ratio 
9/9s. Orig inal model calculations are from Berezhko & 
Volk (|1997|) . 



will only apply for a limited part of the SNR shock sur- 
face. Over other parts the shock is quasi-perpendicular. 
A strongly reduced injection efficiency below that derived 
from injection theory for parallel shocks {rj ^ 10~^ — lO"**) 
along with a commensurate reduction of the particle pro- 
duction will be noticed in such regions (regarding the 
question of injection in general, see Kirk & Dendy 2001 
for a recent review). Correcting for this effect in detail 
is complicated, and will be discussed in a follow-up pa- 
per. However, a rough, but empirically-argued implemen- 
tation of a re-normalisation is achieved by reducing the 
fiuxes F^^ by the ratio 9/Os to a re-normalised flux 

^-13 — 



FZi^Q/Qs, where 9=0.1 is the empirically ex- 



pected value for an average Galactic SNR discussed ear- 
lier (and adopted by DAV) , and 9s is that value predicted 
by the BV model. The re-normalised fluxes are those r- 
subscripted in Table ^, and are rather close to our upper 
limit and the relative CR energy calculated at the present 
age of Tycho. 9r(i) = E„it)/Esn x 9/9s corresponds 
to values of around 10''^ ergs. Fig. ^ graphically compares 
-F2i3 fo'^ the different cases discussed. It appears that cases 
(c) and (d), with a lower injection rate and i?- field, are pre- 
ferred. However, the combinations of B and 77 used here 
are not exhaustive. For example a case assuming higher 
B (perhaps suggested by our inverse Compton interpre- 
tation in section ^ and low injection rate has not been 
tested by BV. More thorough comparisons will be made 



F.A. Aharonian et al.; Tycho's SNR at TeV energies 



later. The uncertainties in scaling parameters, d, Egn and 
n will weaken of course the conclusion favouring cases (c) 
and (d), in which lower injection rates and B fields are 
assumed. The uncertainty in d, where roughly a factor 
of two above the value used here (2.3 kpc) is published 
(Schwarz et al. 1995), will have strong influence. Also, a 
somewhat lower value of n=0.3 cm~'^ derived by Seward 



et al. (1983) for the pre-shock density may also be prefer- 
able to the value given by Smith et al. (1988) used here. At 
this point we would conclude that a general consistency 
with the BV model is obtained, perhaps favouring lower 
injection rates and higher B fields (noting the reduction 
in emissivity between cases (a) and (b), due to a change 
in B field), but note that use of a wider parameter space 
under a dedicated BV calculation is required. 



7. Conclusion 

A search for TeV 7-radiation from Tycho's SNR has been 
performed over two years (1997 & 1998) with the HEGRA 
CT-system. We find no evidence for such emission and the 
upper limit (Scr level) to the TeV flux is estimated at 33 
milli-Crab or 5.78 x 10"^'^ ph cm~^s~'^, a value which is 
a factor ^4 less than that previously published by the 
Whipple collaboration (when assuming a spectral index 
of —1.1 in for the comparison). 

Making use of preliminary RXTE results, it is possible 
to set a lower limit on the magnetic field in Tycho B > 22 
/xG under the assumption that the observed hard or non- 
thermal tail in the X-ray spectrum is attributed to syn- 
chrotron radiation from multi-TeV electrons. A more con- 
servative B field lower limit of ~6 ^G is obtained if we use 
the upper limit to the non-thermal X-ray flux estimated 
by Hwang et al. ( |1998D . It may therefore be premature to 
draw conclusions about the minimum magnetic field until 
full analysis of RXTE data are complete. Interestingly the 
upper limit to the non-thermal energetics of Hwang et al. 
( |1998D is well below that implied by a purely non-thermal 
interpretation of the RXTE and GINGA results. 

Comparisons to a long-adopted model of TeV emis- 
sion fro m the decay of tt° were made (Drury, Aharonian 
& Volk 1994 , denoted DAV), in the initial case assuming 
that Tycho's SNR is well in to the Sedov phase of evo- 
lution. Our upper limit is close to a conservative predic- 
tion of DAV. Uncertainty in the global evolutionary state 
(Sedov, or pre Sedov) of Tycho presents a complication 
upon comparison to the DAV 7r° decay 7-ray predictions. 
As a function of time, advanced kinetic theory (Berezhko 
& Volk 1997, denoted BV) indicates a relatively early 
rise in the 7-ray emissivity of the hadronic channel for 
pre-Sedov epochs, in contrast to the less realistic case of 
single-velocity ejecta, as adopted in the DAV model. We 
find that after adjusting the predictions of BV, firstly scal- 
ing for Tycho's parameters and secondly, re-normalising to 
account for quasi-perpendicular shock directions expected 
in an SNR, a general consistency with our upper limit is 
found for a range of injection rates and B fields. 



We conclude from this non-detection and the interpre- 
tation assuming Tycho to be a source of multi-TeV elec- 
trons and hadrons that Tycho's SNR is yet to be ruled out 
as an average accelerator of Galactic CRs. Our upper limit 
supports the notion that Tycho is in a pre-Sedov evolu- 
tionary state. A more complete investigation will be left 
to a companion paper in which a dedicated calculation in 
the framework of the BV model will be performed for the 
parameters (and their uncertainties) of Tycho's SNR. 

Resulting from a type la supernova, Tycho's SNR in 
principle represents the simplest category in terms of a 
theoretical understanding of the SNR dynamics and inter- 
action with the ISM. Yet, given the complexity in model- 
ing non-thermal processes, the basic consistency between 
experiment and theoretical predictions certainly encour- 
ages future observations of Tycho's SNR and other exam- 
ples of it's class. Because of the \/t dependence on instru- 
ment sensitivity, it is likely that further observations at 
TeV energies of Tycho will have to wait for the next gen- 
eration of telescopes operating in the northern hemisphere 
(e.g MAGIC, Lorenz et al. |1999| ; VERITAS, Krcnnrich et 



al. 1999). With roughly one order of magnitude improve- 



ment in sensitivity for the next generation instruments 
over that currently available, and the fact that current the- 
ory does not leave much room for non-detection at higher 
instrument sensitivity, a decisive test of the question of 
whether or not Tycho's SNR contributes to the Galactic 
CR population at an average level is expected within the 
next few years. 
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